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Using ab initio molecular dynamics we investigate the electronic and lattice structure of AFe2As2 
{A—Ca,, St, Ba) under pressure. We find that the structural phase transition (orthorhombic to 
tetragonal symmetry) is always accompanied by a magnetic phase transition in all the compounds, 
while the nature of the transitions is different for the three systems. Our calculations explain the 
origin of the existence of a collapsed tetragonal phase in CaFe2As2 and its absence in BaFe2As2. 
We argue that changes of the Fermi surface nesting features dominate the phase transitions under 
pressure rather than spin frustration or a Kondo scenario. The consequences for superconductivity 
are discussed. 

PACS numbers: 71.15.Pd,74.62.Fj,61.50.Ks,74.70.-b 



The discovery of iron pnictide superconductors^ with 
critical temperatures Tc up to 57.4 upon doping has 
strongly revived the interest in high-Tj, superconductiv- 
ity. The undoped Fe-based parent compound under- 
goes at low temperatures a structural transition from 
tetragonal to orthorhombic symmetry accompanied by 
a magnetic phase transition to a stripe-type spin-density 
wave state^iii^i^il. While the nature of these two tran- 
sitions is different between LaFeAsO (1111 compound) 
and AFe2As2 (122 compound) with A—{Ba., Sr, Ca), su- 
perconductivity appears in both material classes only 
when the lattice distortion and magnetic ordering are 
suppressed, indicating a strong competition between the 
structural distortion, magnetic ordering and supercon- 
ductivity in iron pnictides. 

Recently, superconductivity in the parent compounds 
1111 and 122 was reported under application of pressure 
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resistivity measurements show superconductivity at « 
12 GPa with Tc = 21 K. In BaFe2As2 superconductiv- 
ity is found to appear gradually with increasing pressure 
while in SrFe2As2 the onset of superconductivity occurs 
abruptly^^. In CaFe2As 2^^i^^'^'* , detailed neutron and X- 
ray diffraction analysis shows that the system undergoes 
a first order phase transition from a magnetic orthorhom- 
bic to a nonmagnetic 'collapsed' tetragonal phase under 
pressure. The possible appearance of superconductivity 
in this 'collapsed' tetragonal phase is presently under de- 
batei^ii^. While various experiments give different values 
of critical pressure9iii^ii^i2ai^i2^i2^i2ii2^ due to the fact 
that the phase transition is sensitive to possible nonhy- 
drostatic pressure effects, Sn content in some samples, or 
the use of single crystals or polycrystalline material for 
structure determination, it is claimed that BaFe2As2 and 
SrFe2As2 do not manifest a 'collapsed' tetragonal phase 
at elevated pressurei^iiiii^i^^iSi. Therefore, the fact that 
structurally similar compounds exhibit phase transitions 
of different nature urgently calls for a theoretical under- 
standing. Moreover, it is still under intensive debate 
which is the driving mechanism of the coUinear stripe- 



type antiferromagnetic ordering; whether the Fermi sur- 
face nesting or the competition of exchange antiferromag- 
netic interactions between the nearest neighbor and next- 
nearest neighbor irons^SiSl^^^^^^^^^^. 

Theoretical work on optimization of cell parameters 
and atomic positions under pressure within the frame- 
work of density functional theory (DFT) has been done 
on ylFe2As2 (A^Ca, Ba, Sr)^i^>ia and LaOFeAs^ML. 
However, by considering Vanderbilt-type ultrasoft pseu- 
dopotentials, Yildirim"^^ obtained a smooth structural 
transition for CaFe2As2 under pressure, observing nei- 
ther a sudden sizable increase of the cell parameters a 
and b nor a strong decrease of the cell parameter c which 
is inconsistent with experimental results^^. Xie et al^ 
optimized within the full potential linearized augmented 
plane wave method (FPLAPW) the orthorhombic lattice 
structure for BaFe2As2 under pressure by relaxing the in- 
ternal parameter zas and the c/a ratio while keeping the 
b/a ratio fixed. This procedure doesn't allow for the de- 
tection of the structural and magnetic phase transitions. 
Opahle et al^ investigated LaOFeAs under pressure and 
found that the system is close to a magnetic instabil- 
ity. Obviously, a complete and unambiguous theoretical 
description of the pressure-induced phase transitions in 
ylFe2As2 (^=Ca, Sr, Ba) is still missing. 

In this paper, we employ the Car-Parrinellc— 
projector-augmented wave'*'^ molecular dynamics method 
at constant pressure^ in order to investigate the 
pressure-induced phase transitions for AFe2As2 (A=Ca, 
Sr, Ba). Since in such transitions the interplay among 
electronic, magnetic and lattice dynamics is essential, a 
combined ab initio DFT with molecular dynamics ap- 
proach as used in the present work is very suitable. In 
such a procedure, at each pressure value a full, unbi- 
ased relaxation of all lattice and electronic degrees of 
freedom is performed in time steps of 0.12 fs at zero 
temperature. We used 4x4x4 fc-points in doubled 
(\/2 X \/2 X 1) unit cells. We used high energy cutoffs of 
612 eV and 2448 eV for the wave functions and charge 
density expansion respectively. The total energy was con- 
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FIG. 1: (Color online) Calculated changes of (a) volume, c/a 
ratio and magnetization (inset), (b) lattice parameters as a 
function of applied external pressure normalized to the critical 
pressure for CaFe2 As2 . The phase boundary between F rarara 
and I A/mmm is indicated by the vertical dashed line. 
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FIG. 2: (Color online) Calculated changes of (a) As-Fe-As 
angles, and (b) Fe-Fe, Fe-As bond length as a function of 
applied external pressure normalized to the critical pressure 
for CaFe2As2. The phase boundary between Fmrnm and 
I A/mmm is indicated by the vertical dashed line. 



verged to less than 0.01 meV/atom and the cell parame- 
ters to less than 0.0005 A. The 3s3p3d (4s4p4d/5s5p5d) 
states in Ca (Sr/Ba) and the 3d4s4p states in Fe and As 
are treated as valence states. We checked our calcula- 
tions with the FPLAPW method as implemented in the 
WIEN2k code^. Very good agreement is found between 
these two methods. The Perdew-Burke-Ernzerhof gener- 
alized gradient approximation (GGA) to DFT has been 
used. 

Our findings can be summarized as follows: We show 
that CaFe2As2 undergoes a first order phase transition to 
a collapsed tetragonal phase with relative changes of lat- 
tice parameters, bond lengths and angles agreeing very 
well with the experimental observations. In particular 
our calculations can account for the measured expan- 
sion along the ab plane^^ at the critical pressure, which 
was not obtained in previous calculations. We also find 
an abrupt disappearance of magnetization at the criti- 
cal pressure. For SrFe2As2 and BaFe2As2, where less is 
known about the details of the lattice structure and mag- 
netization changes under pressure, we obtain a simulta- 
neous structural (orthorhombic to tetragonal) and mag- 
netic phase transition at high pressure. However, we ob- 
serve a weak first-order phase transition for SrFe2As2 and 



a continuous phase transition for BaFe2As2 in contrast 
to the strongly first order phase transition in CaFe2As2. 
We explain the microscopic origin of these differences in 
terms of the different Fermi surface behavior under pres- 
sure. Finally, we argue that the existence of a magnetic 
transition from a striped AF state to a non-magnetic 
state under pressure is mainly driven by changes of the 
Fermi surface nesting rather than spin frustration or a 
Kondo scenario'*^. 

In Figs. [1] (a) and (b) we present the calculated changes 
of the volume and lattice parameters of CaFe2As2 as a 
function of pressure. The volume and the lattice param- 
eter c decrease gradually with increasing pressure and 
show a discontinuous shrinkage at the critical pressure, 
where the system undergoes a structural phase transi- 
tion from orthorhombic symmetry to a volume 'collapsed' 
tetragonal symmetry. Our results are in very good agree- 
ment with experimental data^^ with a volume collapse 
of Al/'h « 4.1%, AF^'^P « 4%. Moreover, the mag- 
netization sharply goes to zero as shown in the inset of 
Fig. [l](a). Surprisingly, however, the lattice parameters 
a and b show abrupt expansions at the phase boundary 
although otherwise they show a monotonous compres- 



3 



dxz/dyz 




K 



Z r X P NO DOS 15 




DOS 



FIG. 3: (Color online) Comparison of the CaFe2As2 band 
structure and DOS in the 'collapsed' tetragonal phase ob- 
tained from (a) our optimized lattice structure at P/ Pc = 1.04 
and (b) the experimental lattice structurei^ at p'^^p/p^^p ~ 
1.6. In (a) the total DOS for the high T tetragonal phase is 
also shown by the dashed (red) line. 



sion in both phases. Such expansions, as also observed 
in experiments^^, were not detected in previous DFT cal- 
culations and have a fundamental physical origin. In the 
orthorhombic phase below the critical pressure, the sys- 
tem shows collinear AF ordering within ab plane. At 
the critical pressure, the AF ordering is destroyed and 
the Pauli principle becomes more effective due to the 
appearance of more parallel Fe spins giving rise to a non- 
negligible lattice expansion along a and b^. As a con- 
sequence, the Fe-Fe bond length increases at the phase 
transition as observed in Fig. [2] (b). In contrast, the As 
atoms move towards the Fe plane at the phase transi- 
tion and the Fe-As bond length suddenly decreases with 
a relative shrinkage of Ad^p^_j^^ w 1.3% in very good 
agreement with Ad'^p^^^^ ~ 1.2%. Accordingly, the an- 
gles of the tetrahedra are suddenly shifted away from the 
ideal tetrahedral value of 109.47° at the phase transition, 
while in the orthorhombic and in the high temperature 
tetragonal phase they are much closer to the ideal value 
(Fig. [5] (a)). The strong distortion of the tetrahedron 
enhances the crystal field splitting between (i^2_j,2, d^z, 
dyz and d^^, dxy (orbitals given in the a; || a, y || 6, z || c 
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FIG. 4: (Color online) Calculated changes of (a) volume, c/a 
ratio and magnetization (inset), (b) lattice parameters and 
(c) As-Fe-As angles as a function of applied external pressure 
normalized to the critical pressure for SrFe2As2. The phase 
boundary between Fmmm and 1 4/mmm is indicated by the 
vertical dashed line. 



local reference frame). This can be seen in the density 
of states (DOS) shown in Fig. [3] (a) where the two peaks 
just below (of dj;2_y2^ d^z, dyz character) and above (of 
dz2, dxy character) the Fermi level in the high temper- 
ature ambient pressure tetragonal phase (dashed (red) 
line) are further separated in the high pressure volume 
'collapsed' tetragonal phase (solid (black) line) . 

One aspect of the calculations to be considered is that 
the calculated critical pressure Pc = 5.25 ± 0.25 GPa is 
an order of magnitude larger than the experimental crit- 
ical pressure P^'^p « 0.3 GPa. The overestimation of 
Pc in our calculations is a consequence of the well-known 
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FIG. 5: (Color online) Calculated changes of (a) volume, c/a 
ratio and magnetization (inset), (b) lattice parameters and 
(c) As-Fe-As angles as a function of applied external pressure 
normalized to the critical pressure for BaFe2As2. The phase 
boundary between Fmmm and 1 4/mmm is indicated by the 
vertical dashed line. 



overestimation of volume and magnetic moment by GGA 
in the Fe pnictides at ambient pressure^i^. Therefore 
in order to reach the experimental structure conditions 
we have to go to higher simulation pressures. This pro- 
cedure is proven to be valid since, as discussed above, 
we can reproduce the experimentally reported structural 
features of the phase transition in CaFe2As2, as well as 
the electronic properties in both orthorhombic phase"^? 
(not shown) and 'collapsed' tetragonal phase as shown 
in Figs. [3] (a) and (b) where we present the comparison 
between the electronic structure obtained from the Car- 
Parrinello-derived lattice structure (at P/Pc = 1.04) and 




FIG. 6: (Color online) Fermi surface cuts centered at (0,0,fcz) 
with kz =0, tt/c and 27r/c and \kx\, \ky \ < n/a at three pres- 
sure values for CaFe2As2. The Fermi surface sheets are shown 
in black (blue). The grey (orange) sheets result from translat- 
ing the sheets centered at (tt/o, 7r/a, fe^) by (— yr/a, — 7r/a, 0) 
in order to see the nesting effects at (0,0,fcz). Note, that the 
low P Fermi surface should not be compared to ambient P 
high T Fermi surface as we always work at zero temperature. 
Here the pressures are normalized with respect to the critical 
pressure Pc for CaFe2As2. 



that obtained from the experimentally measured one (at 
0.47 GPa, pexp^pexp _ observe that the shapes 

of the total DOS are almost the same and only slight 
differences are found in details of the band structure. 
Therefore we can conclude that, although our calcula- 
tions always overestimate the critical pressure, the phase 
transition and the physical properties in both phases can 
be quantitatively captured by our analysis. 

In Fig. m and Fig. O we show the results for SrFe2As2 
and BaFe2As2. We find that the structure and the mag- 
netic phase transitions under pressure still occur simul- 
taneously in these two compounds, but the nature of the 
phases is distinctly different from the CaFe2As2 case. 
SrFe2As2 shows smaller magnetization, volume, lattice 
and angles jumps at the critical pressure (see Fig. |4]) 
compared to CaFe2As2 while we can hardly detect any 
discontinuity in BaFe2As2 (see Fig. [SJ. The simultane- 
ous and abrupt change of the lattice structure and the 
magnetization in SrFe2As2 is again consistent with recent 
experimental findingsilii^ while our results in BaFe2As2 
have not yet been completely confirmed experimentally. 
Here, the pressure is again normalized to the critical sim- 
ulated pressure P^ = 9.5 ± 0.5 GPa for SrFe2As2 and 
P^' = 17.5 ± 0.5 GPa for BaFe2As2. The calculated crit- 
ical pressure for SrFe2As2 and BaFe2As2 is about 2-3 
larger than the experimental one, while for CaFe2As2 it 



FIG. 7: (Color online) Fermi surface cuts centered at (0,0,^2) 
with kz = 0, tt/c and 2n/c and \kx\, \ky\ < n/a at three 
pressure values for BaFe2As2. The color coding is the same 
as in Fig. |6] Here the pressures are normalized with respect 
to the critical pressure P" for BaFe2As2. 
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FIG. 8; (Color online) Total DOS for CaFe2As2 (continuous 
line) and BaFe2As2 (dashed line) at the critical pressure re- 
gion. 



is around 20 times larger. This can be understood from 
the already mentioned overestimation of magnetic inter- 
actions and volumes in GGA calculations. Due to this, 
one would expect an (approximately constant) shift be- 
tween the calculated and measured Pc, thus resulting in 
the large deviation for CaFe2As2 (where Pc is small) and 
the better agreement for SrFe2As2 and BaFe2As2. 

In order to understand the microscopic origin of the 
differences in behaviour between SrFe2As2, BaFe2As2 
and CaFe2As2, we performed non-spin-polarized calcu- 
lations with the optimized lattice structures and ana- 
lyzed the band structures and Fermi surfaces. Below the 
critical pressure, we considered an average of lattice pa- 
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FIG. 9: (Color online) Calculated noninteracting suscepti- 
bility Xo{<l) at g = (tt, 7r,0) within the constant matrix el- 
ement approximation as a function of the normalized pres- 
sure with respect to the critical one (in arbitrary units). This 
magnitude measures the response for the observed stripe-type 
antiferromagnetic ordering. The three panels correspond to: 
(a) CaFe2As2, (b) SrFe2As2, and (c) BaFe2As2. The phase 
boundary between Fmmm and I i/mmm symmetries is in- 
dicated by the vertical dashed line. 



rameters a and 6 as a single lattice parameter in order 
to make the comparison of band structures of different 
phases (orthorhombic and tetragonal) possible within the 
same Brillouin zone. It has already been noticed that the 
tiny differences between a and b has no appreciable in- 
fluence on the electronic structure^. 

In Figs. [6] and [7] we show kz =0, tt/c, and 27r/c Fermi 
surface cuts of CaFe2As2 and BaFe2As2 at three different 
pressures. As a function of pressure, the changes of the 
Fermi surface in both compounds is remarkably differ- 
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ent. Below Pc, CaFe2As2 shows with increasing pressure 
gradual shrinkages of the nested inner hole and electron 
cylinders around (0,0,^^) and (tt/o, 7r/a, /c^) respectively 
leading to a reduction of the nesting features at each kz 
(compare in Fig.[6]thc rows at P/Pc=0.1 and P/Pc^O.95) 
and consequently a continuous reduction of the magneti- 
zation in the orthorhombic phase (inset of Fig.[T] (a)). As 
pressure increases to the critical region, the residual steep 
DOS at the Fermi level (Fig. [8]) still forces the system to 
find a way to alleviate such an instability. However, the 
nested inner cylinders become too small to remove it by 
magnetic ordering. A sudden distortion of the tetrahe- 
dron is realized to lower the electronic energy by further 
splitting the Fe 3d orbitals (see the dashed (red) DOS hue 
in Fig. [3] (a)) at the price of increasing the lattice energy. 
Thus, as shown in Fig. [SI P/Pc=l-05 row, right above the 
critical pressure all the cylinders around (0,0,^^) com- 
pletely vanish due to the splitting of Fe 3d orbitals by 
lattice distortions in contrast to BaFe2As2 (see below), 
and accordingly, the magnetic ordering abruptly disap- 
pears (inset of Fig. [T] (a)). 

In contrast, the BaFe2As2 Fermi surface shows a 
shrinkage of the hole cylinders around (0,0, fc^) into 
hole pockets while the Fermi surface sheets around 
{TT/a,n/a,kz) are changed insignificantly (see Fig. [7]). 
Therefore, as the shrinkage occurs under pressure, the 
vanishing of Fermi surface nesting at the lower ft^-plane 
is partially compensated by the appearance of new nest- 
ings at a higher fc^-plane (e.g. Fig. [7] at P/P^'=0.17 
and P/P"=0.97), which prohibits the sudden change of 
magnetization since the nesting remains significant be- 
low the critical pressure (e.g. Fig. [7] at kz = 27r/c and 
P/P"=0.97). Since lowering the total energy by mag- 
netic ordering always prevails over abrupt lattice distor- 
tions, the sudden jump of the angles is suppressed as 
shown in Fig. [5] (c). The slow change of the Fermi sur- 
face topology as shown in Fig. [7] is consistent with the 
continuous reduction of magnetization as observed in the 
inset of Fig. [H] (a). 

Fig. [9] presents the calculated noninteracting suscep- 
tibility xo{q) at 5 = (tt, 7r,0) as a function of pres- 
sure within the constant matrix element approxima- 
tioni^i2^i2^i2^ for all the three 122 compounds. The re- 
lation between the change of the Fermi surface nesting 
features at g = (tt, tt, 0) and the phase transitions un- 
der pressure is clearly quantified with this magnitude. 
Xo{q) at g = (tt, TT, 0) is the response for the observed 
Fermi-surface-nesting-induced stripe-type antiferromag- 



netic ordering and while it shows a continuous decrease 
at elevated pressures in BaFe2As2, abrupt reductions at 
the phase transition are detected in both CaFe2As2 and 
SrFe2As2, indicating that the phase transition for these 
two systems are of first order where a sudden disappear- 
ance of the Fermi surface nesting occurs. Moreover, we 
observe that the volume and lattice parameter discon- 
tinuities (Fig. m (a) (b)) and the tetrahedron distortion 
(Fig.[4](c)) right above the phase transition for SrFe2As2 
are of the same nature as in CaFe2As2 but less pro- 
nounced. We conclude that the phase transition is of 
weak first order in SrFe2As2. 

Finally, we find that the ground state energy of the 
stripe-ordered AF state is always far below that of 
the checkerboard-ordered one in the orthorhombic phase 
while that of the high-pressure tetragonal phase is non- 
magnetic, indicating that the phase transition under 
pressure for the three compounds is not driven by the 
competition between different magnetic ordering, i.e., 
spin frustration. Furthermore, the tiny crystal field split- 
tings of the Fe 3d states and strong itinerancy of all Fe 
3d electrons does not support the Kondo scenario^ for 
the phase transitions under pressure. 

In conclusion, we investigated the phase transitions in 
the iron pnictide 122 compounds under pressure within 
the framework of ab initio molecular dynamics. The re- 
sults on CaFe2As2 agree very well with the experimental 
findings. We show that magnetic and structural phase 
transitions also appear concomitantly in SrFe2As2 and 
BaFe2As2. While they are weakly first order in SrFe2As2, 
the phase transition is continuous in BaFe2As2 which 
may be closely related to the different behavior of enter- 
ing superconductivity between SrFe2As2 and BaFe2As2 
observed experimentally^^. The remaining hole pocket as 
well as nesting features in BaFe2As2 may be related to 
the highest transition temperatureSS observed among 122 
compounds due to the remaining spin fluctuations. Since 
the topology of the Fermi surface in the high pressure 
phase is quite different among these 122 compounds and 
LaOFeAs^, the pairing mechanism may differ among 
these compounds although they belong to the same fam- 
ily of superconductors. 
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